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Abstract 

Transition from quantum to semiclassical behaviour and loss of quantum co- 
herence for inhomogeneous perturbations generated from a non-vacuum initial 
state in the early Universe is considered in the Heisenberg and the Schrodinger 
representations, as well as using the Wigner function. We show explicitly that 
these three approaches lead to the same prediction in the limit of large squeez- 
ing (i.e. when the squeezing parameter |rfc| — ► oo): each two- modes quantum 
state (k, — k) of these perturbations is equivalent to a classical perturbation 
that has a stochastic amplitude, obeying a non-gaussian statistics which de- 
pends on the initial state, and that belongs to the quasi-isotropic mode (i.e. 
it possesses a fixed phase). The Wigner function is not everywhere positive 
for any finite r^, hence its interpretation as a classical distribution function in 
phase space is impossible without some coarse graining procedure. However, 
this does not affect the transition to semiclassical behaviour since the Wigner 
function becomes concentrated near a classical trajectory in phase space when 
|rfc| — > oo even without coarse graining. Deviations of the statistics of the 
perturbations in real space from a Gaussian one lie below the cosmic variance 
level for the iV-particles initial states with N = iV(|k|) but may be observable 
for other initial states without statistical isotropy or with correlations between 
different k modes. As a way to look for this effect, it is proposed to measure 







the kurtosis of the angular fluctuations of the cosmic microwave background 
temperature. 

PACS Numbers: 04.62. +v, 98.80. Cq, 98.80.Hw 
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1 Introduction 



The undeniable success of the theory of the generation of cosmological perturbations 
from quantum vacuum fluctuations during an inflationary era in the early Universe 
was its ability to correctly predict in advance three main properties of the recently dis- 
covered (non-dipole) angular fluctuations AT ^ of the cosmic microwave background 
(CMB) temperature. These predictions are: a) the multipole power spectrum of the 
fluctuations that is approximately flat for low multipoles (Cj oc j^jy for 2 < 1 < 30); 

b) the right location and amplitude of the first acoustic (or "Doppler" ) peak of 4rr at 
I = 200 — 250; c) the Gaussian statistics of the fluctuations (see, e.g., the 4- year COBE 
results for the spectrum |jj and statistics of the fluctuations with I < 20, and the 
results of the Saskatoon-95 experiment for I = 80 — 350 and the CAT experiment 
for I = 400 — 600 J3J). Of course, this does not exclude alternative theories of the early 
Universe (e.g., topological defects) provided they are able to produce approximately 
the same observational predictions. However, it is clear that we can and even must 
seriously investigate quantum-gravitational effects in the early Universe because of 
the actual observability of their consequences at present time. 

At this point, one of the most important and interesting questions of this theory 
arises. The initial perturbations generated during an inflationary stage were quan- 
tum, relativistic and gravitational since their dimensionless amplitude contains all 
three fundamental constants: the Planck constant H, the light velocity c and the 
Newtonian gravitational constant G (below we will take c = h = 1 unless otherwise 
stated). The fluctuations -J- measured today, on the other hand, are certainly classi- 
cal (though stochastic). Thus, some mechanism describing the quantum-to-classical 
transition of the perturbations is required. What was explicitly or implicitly used 
in the pioneering papers where the spectrum of tensor perturbations (gravitational 
waves) and scalar (adiabatic) perturbations || generated from vacuum fluctu- 
ations was first correctly derived, was a kind of equivalence between the strongly 
squeezed vacuum fluctuations and classical perturbations with a stochastic Gaussian 
amplitude and a fixed temporal phase in the Heisenberg representation (see [pj for 
a detailed discussion and explanation of this equivalence). The same applies to the 
Schrodinger representation, see || where the quantum cosmological approach based 
on the Wheeler-de Witt equation was used, it reduces to the Schrodinger representa- 
tion for each mode after separation of variables and introduction of the quasi-classical 
time. Finally, this equivalence can be seen as well using the Wigner function formal- 
ism |?|, |9|. An essential property of this approach is that this type of transition from 
a quantum to a classical description of the perturbations, including the complete 
loss of quantum coherence, does not require any actual consideration of very small 
interactions between a given mode of perturbations and other degrees of freedom 
("environment") in the Universe. Though these interactions do of course exist and 
do lead to the transformation of an initial pure quantum state into a mixed state, their 
actual structure is unimportant for observable quantities as far as the exponentially 
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small decaying mode of perturbations may be neglected (note however that it must 
be taken into account if one wants to calculate the entropy of cosmological pertur- 
bations flOf). This is precisely why this transition was called "decoherence without 
decoherence" in J7|. Furthermore, though the quantum coherence disappears, not 
all coherence is lost: indeed, the equivalent classical stochastic perturbation remains 
partly correlated (coherent) at the classical level in the sense that it has a fixed tem- 
poral phase corresponding to the fact that it belongs to the special quasi-isotropic 
mode of perturbations. Remarkably, this coherence has observable consequences: it 
results, in particular, in the appearance of multiple peaks (called acoustic, or Doppler 
peaks) in the multipole power spectrum of ^f- with periodic spacing A/ = n ll0 ~^ ec 
in the case of adiabatic perturbations. Here v s is the sound velocity in the relativis- 
ts plasma consisting of coupled baryons and photons just before recombination, r\ 
resp. t] rec are the conformal times today resp. at recombination. A similar effect 
(though of significantly smaller amplitude) takes place for CMB angular temperature 
fluctuations produced by primordial gravitational waves background @ (v s — 1 in 
this case). 

This approach however, especially that one based on the Wigner function for- 
malism, was strongly criticised in [|llj] on the ground that it might be applicable to 
vacuum initial conditions only and not to more general ones. Hence, we adress in 
this paper the problem of the quantum to classical transition of quantum cosmological 
perturbations in the case of non- vacuum initial states. By an "initial state" of a given 
Fourier mode k of a quantum field we mean its state in the WKB regime deep inside 
the Hubble radius at the inflationary stage, long before the Hubble radius crossing, 
i.e. when k ^> aH where k = |k|, a(t) is the scale factor of the Friedmann-Robertson- 
Walker (FRW) isotropic cosmological model, H = a/ a is the Hubble parameter and 
the dot means a time derivative with respect to the cosmological time t. 

Of course, a non-vacuum initial state contradicts the whole spirit of the maximally 
symmetric initial state of the Universe which lies at the heart of the inflationary sce- 
nario. Thus, our "wish" would rather be to show that such an assumption is in 
conflict with observational data. Surprisingly, as we will see below, this turns out to 
be more delicate. Non-vacuum initial states are typically not self-consistent in the 
simplest variants of the inflationary scenario since they lead to a large energy density 
of inflaton field quanta which is not of the type of a cosmological term and it would 
contradict the assumption of inflation at time rj especially if we push t]q to — oo. 
However, they could appear in more complicated, double inflationary models for per- 
turbation modes which cross the Hubble radius three or more times during inflation 
if we additionally assume (without justification) that these modes immediately deco- 
here (e.g., in the particle number basis) once they reenter inside the Hubble radius. 
Though the latter process is hardly physically viable, still it cannot be completely 
excluded. Furthermore, it is an important question of principle to study whether 
the possibility of a classical description of highly squeezed perturbations is generic or 
whether it is restricted to vacuum initial states only. 
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Particle creation from non- vacuum initial states was considered in many papers; 
in the cosmological case - beginning with the pioneering paper by Parker |T2 
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the recent review ]T3[ for many references. However, the primary quantities studied 
in these papers were the number of created particles and the average value of the 
energy-momentum tensor of a quantum field on a FRW background. In contrast, we 
are mainly interested in the field fluctuations themselves since these are the quantities 
(and not the number of particles) that are directly observable in the case of cosmolog- 
ical perturbations. This refers both to scalar metric, or density, perturbations and to 
tensor metric perturbations or gravitational waves. Naturally, the resulting quantum- 
to-classical equivalence that, as we argued in 0, is the reason for the possibility of a 
classical description of the perturbations in the present Universe takes place for these 
quantites, too. 

Thus, the results of paper [0 are generalized to the case of non-vacuum initial 
conditions. In order not to introduce any preferred direction in space, we will consider 
only reflection symmetric (r — > — r, k — > — k) quantum states. We show that a 
wide class of initial states become quasiclassical in the limit of high squeezing. This 
is shown in section 2 using the Schrodinger representation, and in section 3 using 
the Heisenberg representation. In section 4 we derive conditions necessary for the 
semiclassical behaviour of quantum scalar (adiabatic) perturbations and gravitational 
waves on a FRW background when expressed directly in terms of rms values of the 
perturbations. The Wigner function approach to this problem is studied in section 

5. The specific, but important case of a thermal initial state is considered in section 

6. Finally, section 7 contains conclusions and a discussion of a possible observational 
confirmation or rejection of a non-vacuum initial state. 



2 Schrodinger representation 

In this section, we are going to compute the wave function of different non-vacuum 
states in the Schrodinger representation. Let us first review the basics about the 
dynamics of the system that we consider here, namely that of a real massless scalar 
field on a flat FRW universe. Actually, the field can be massive and all the arguments 
presented here will go through. The space-time metric has the form 

ds 2 = dt 2 — a 2 {t)5ijdx % dx^ , i,j = 1,2,3. (1) 

Let us remind how the dynamics of this system will lead to the appearance of 
squeezed states. We first write down the Hamiltonian H in terms of the field y = a<p 
and the conformal time rj = f ^y. The following result is then obtained 

H = J d 3 xH(y,p,d iy ,t) 

= ~ / d 3 %(k)pt( k ) + **y(k)yt(k) + - (y(k)pt( k ) + p(k)j/t(k))] (2) 
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where 

dC(y,y') , a' 

P = — tt-, — = y — y 3) 
ay' a 

and a prime stands for derivation with respect to conformal time. Here the following 
Fourier transform convention is used: $(k) = (27r)~ 3 / 2 / $(r)e~ lkr (i 3 r for functions 
as well as for operators. Due to reality of the field y, we have that y(k) = y*(—k), 
resp. y^{— k) for operators. Therefore, any classical field configuration is completely 
specified by giving the Fourier transforms in half Fourier space. Though generally this 
doesn't have to be the case when the fields are quantized, we will deal here with states 
which are invariant under the reflection k — > — k as mentioned above. The system can 
be described with the help of the complex valued field modes fk(f]) with $lf k = f k i 
and ^sfk = fk2, fk(vo) = l/V^k (we adopt a similar notation for all quantities), and 
momentum modes g k (i]), Qkijlo) = \/f > typically used in the Heisenberg representation 

y(k,rj) = f k ( V ) a(k,r] ) + f*( V ) a\-k,r] ) , (4) 
= -i (g k (v) o(k,?7o) - g* k {rj) ^(-k,^)) . (5) 

In terms of f k and g k , the important quantities 7^ and F{k) are given by the following 
expressions 

lk = ~F = 2L^F (1 ~ l2F{k)) ' F{k) = fkl9k2 ~ fk29kl ' (6) 

The functions f k satisfy the equation 

ti+(k 2 -^jfk = ^ (7) 

and 

9k = i(fL- ~fk) ■ (8) 



The couple of functions f k , g k represent only three independent real functions 
as they satisfy a Wronskian constraint. The dynamics of the system can also be 
conveniently parametrized by the following three functions: the rotation angle 9 k , the 
squeeze angle (f> k and the squeeze parameter r k . When \r k \ — > 00, the system is said 
to be highly squeezed. In that case also \F k \ ^00. As a result of the time evolution 
of our system, we have 

S a f (±k,r7o) = 9k(v) y{TKvo) ~i fk(v) p{tKVo) , 

S a(±k, Vo ) S' 1 = ^)y(±k,i to )+i/;(i 7 )p(±k, % ) (9) 

where S is the unitary evolution operator from time 770 to time 77. When working in the 
Schrodinger (complex) coordinate representation, we will often adopt the following 
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compact notation y(k, r) ) = y , p(k,r) ) = Po = _i ^> y(~KVo) = Vo, p(-KVo) = 
Pq = —i-^ where yo and y$ are considered as independent variables. 

As a result of the coupling with the gravitational field, which plays here the role 
of an external classical field, the state |0, 77)5 which is the vacuum state of the field 
at some given initial time rj , a(k, rj ) |0, rj Q ) = Vk, will no longer be the vacuum 
at later times rj. Indeed field quanta are produced in pairs with opposite momenta, 
we get a two-modes squeezed state. Though the state at late time is no longer 
annihilated by the annihilation operators a(k), it is still annihilated by the following 
time-dependent operator 

{y(k, no) + P(k, r/o)}|0, v)s = 0. (10) 

In the infinite volume case, equation flTU| ) translates into the following functional 
differential equation in the Schrodinger (complex) coordinate representation 

{y^)+%\v) I ^^ ) }\0,v)s = 0. (11) 

Therefore the state at late time is still Gaussian, its wave functional in the infinite 
volume case (the quantum variables are the functions y(k) of the continuous variable 
k) being given by 

V [y(k, T) ),y(-k, rj )} = Af k exp (-- J d 3 k 7fc y(k, r) ) y(-k, r) )\ . (12) 

For our purposes however, it is enough to consider our system enclosed in a finite 
volume and we then get as wave function of our system an (infinite) product of 
functions 

^-(-M {1 _ !2W} ) (13) 

for each pair k, — k. The commutators and correlation functions of Fourier transform 
operators are then c-number functions and not distributions. The time dependence 
of \1/ is through f^, F(k), and iV^. We note the persisting symmetry under reflections 
in k-space, it is present in the initial state and the evolution doesn't spoil it. 

As a first generalization of the vacuum state at time r] , we consider a two-modes 
coherent state \a(3, rj )s satisfying 

a(k,r] ) \af3, T] ) s = -y(k,r] ) + -^=p(k,7] )j \a/3, 7] ) s 

= a I a/3, 7] )s = {y - (y)o + (p)o\ \a/3, 770)5 (14) 
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where (y) and (y) are the average values of y(k, r/ ) and p(k, 77 ) at time rj , and it 
satisfies also 

a(-KVo) \otP, Vo)s = {J |?/(-k,?7o) + -^L=p(-k, 77 )} |a/3, 770)5 



= (3 \a/3, Vo ) s = {^(yy, + -J=( P y } 770)5 . (15) 

It is clear that the time-evolved state \aP, rj) s satisfies the following equations at any 
later time 77 

S a(k,77o) S 1 - 1 \a(3, r}) s = a \a/3, rj) s = ( g* k (rj) (y) +i ft(rj) (p) ) |a/3, r)) s 
S a(-k, 770) V)s = P \<xP, 77)5 = ( ^(77) (7/)* + i f* k (v) ( P )* ) |a/3, 

where (y), (p) are average values at time 77. It follows from ( pl4| , |T6| ) that (y) is the 
solution of eq. (^) with the initial conditions 7/(770) = (y) , p(i] ) = (p) , and 
(p) = (y)' ~ ~(v)- Eqs. flTfi|) are easily solved in the Schrodinger representation 
resulting in the following wave function ^ap(yo, Vo, v) °f the two-modes coherent state: 

*a/3(yo, 77) = MMo),2/(-k, 77 )|a/?, 77)5 

= iV fc exp[ - lk (y - (t/))(t/* - (y)*) + i(y* (p) + y (p)*) ] (17) 
= iV fc exp[ - 7fc \y - (y)\ 2 + i( y* (p) + t/ (p)* ) ] ■ (18) 

We could equally well arrive at that result working in the real and imaginary co- 
ordinate representation. Indeed, we can separate the variables y and p in real and 
imaginary part and consider the corresponding operators in the Schrodinger repre- 
sentation. Defining S a(±k, 770) S* -1 = (21(77) ±20,2(77), the above equations can be 
written as : 

01(77) \ap, v)s = {g* k (v) (Vi) +ifk(v) (Vi)} \a/3, v)s 

1 d 

01(77)^(7/0,7/0,77) = gl[T]) 7/01 + fkW^Q-^ ^*p(yo,ya,V) , (19) 

02(77) \ap, 77)5 = {^(77) (t/ 2 ) + ifk(v) (2/2)} \aP, 77)5 

1 9 

o 2 (?7) *a0{Vo,yo,r)) = 9k(V) fe + Ml)^ ~ *ap(yo>Vo>V) • ( 20 ) 

Solving this differential system, we find that the wave function can be written as: 
^apivOiVoiV) = 0(2/oi, 77)0(7/02,77), where : 

0(2/oi, ^) = \/iVfcexp[ -7^ (7/01 - (t/i)) 2 + 72 7/01 (pi) ] , (21) 
in complete accordance with ( |T7j|i"8| ). 
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We are further interested in states which, at some initial time r] , contain some 
arbitrary number iV of particles with momenta k and — k and we denote these states 
by |iV k ,AL k , T] ) s = \N, Vo)s- 

\N, V o)s = ^at iV (k,r ?0 )a^(-k,r ?0 ) |0, Vo ) s . (22) 

We deal with modes initially deep inside the Hubble radius (k ^> a' (rj ) / a(rj )) and use 
the WKB solution for these modes. The time-evolved state defined in (1221) becomes : 



\N, V ) s = ^(S a\k, m ) S-'f (S a^-k,^) S-'f \0, rj) s . (23) 

In the Schrodinger coordinate representation, the following differential equation relat- 
ing the N-particles two-modes wave function to the vacuum two-modes wave function 
is obtained: 

1 / d \ N f d \ N 

^n(vo, y* , v) = j^fk N (it vo - -^A hi vl - ^ J Mvo, vl, v) (24) 

with 

^o(yo,yo,v) = (V^IAl) 1 ex P ("Tit yo Vo) 

= N k exp (-7 fc y y* Q ) . (25) 

The action on \l/ of the operator inside the right parenthesis can be easily derived 
using the relation 7^ + 7^ = \fk\~ 2 - Indeed: 

* N = Ti(fi) N *(rty°-i^) yf eM-ik yo y* ) • (26) 



We make further a binomial expansion of the first parenthesis : 

^ = M (|) N k Y^CUll yo) p ( ) {^exp (- 7fe y y*o)} • (27) 



Ik/ p=0 

It is now convenient to use the following trick, namely to replace (7^ yo) p by the 
equivalent expression: 

( 7fc * yo) p = exp( 7 ^ yo y* ) (~) {exp(- T £ V° Vo)} ( 2 §) 



As a result we get 

' f , \ N N ( r> \ P 

{exp(-7fc yo y* )} x 



d \ N ~ P 



{exp(-7 fc y y*)}, (29) 
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where we recognize the Nth derivation of a product of functions. The Leibnitz formula 
gives: 



A? 



J_(fk 

N\ [ft 



N 



iV fc exp(7fc y y* ) 







dyo. 



N 



N 



exp 



l/o 
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I/*IV V^o 



*7V 

2/o ex P 



1 2/o I 



I /* 



(30) 



The expression inside the square brackets is just a Laguerre polynomial as can be 
checked using Rodrigue's formula: 



L N (w) 



(31) 



Taking w(y*) = — jj^Tp 2/o, where ?/o is constant, it follows that the expression inside 
the square brackets is just Ln The ^ na ^ res ult is therefore: 



TV 



*tv(j/o, 2/0, ?7) = f j ^V J #0(2/0, 2/o, ^ • 



(32) 



This calculation can also be done differently, in a way which takes advantage of basic 
results pertaining to the familiar one-dimensional harmonic oscillator. Indeed, let us 
separate the complex variables and their corresponding operators in real and imag- 
inary parts, as we did already previously. The definition of the state then becomes 



\N, v)i 



-^ (S at (k, 7ft) S- 1 S a\-k, Vo ) S'^Mi 
— (of + 4 2 ) JV |0,77) S 



^(EC4 2M 4 2iV - 2M )|o^) s . 

iV - \M=0 ) 



In this representation, the differential equation 
^(2/0,2/0,77) 



becomes: 



2N-2M 



yrtv) e cjf U 2/01 - 

( 7 * ym - si 

TV / 

m=o V 



1 d 



2M 



2dy 01 



N, 
M 



#0(2/0, Vo,v) 

o \ 2M 
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2 dy i . 

2N-2M 



| e -7fe2/giJ x 



7 fe 2/02 - - 



-7fc 



3/02 \ 



2 9?/02 



(33) 



(34) 
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Knowing the harmonic oscillator eigenfunctions, it is straightforward to find: 



N 




pM TT 

'-'AT a 2M 




2N-2M 



V02 
\fk\ 



(35) 



It is now possible to show that this expression is equal to (B^). Taking into account 
the relation between Hermite and Laguerre polynomials, we write: 



N 



< / p \ f JV i/2 



im 2 . 



X 



- x - ]L > 2X - 2M (N-M)\L ( -p M ( 2/02 




AT! I /• I 1 ^ 



(36) 



where 8.974.4, p. 1038 of [|14]] has been used in the last step of the above calculation. 



The two expressions (|32| , |35D are two equivalent ways to write down the wave 
function in the coordinate Schrodinger representation. They represent, together with 
(|I~7D, the basic result of this section. We note that the wave function (|32|) is a function 
of |y(k)| and therefore, like the vacuum wave function, it is invariant under reflection 
k — > — k. As a result, it is clear that all correlation with an odd number of points will 
vanish, the distribution has no skewness. On the other hand, it is of course no more 
a Gaussian, a fact which will be important regarding the positivity of the Wigner 
function. The transition to semi-classicality can now be seen in this representation. 
The condition for a semiclassical behaviour is that the phase of the wave function 
should be large for a typical value of \y\ 2 ~ (\y\ 2 ) = (2iV + l)|/fc| 2 . This results in 
the inequality (2iV + l)\F k (r])\ ^> % (here we restored % for clarity). Note that, as one 
could expect, the latter condition is more easily achieved for large N. The probability 
distribution for 7/0 a t time 77 is given by: 



pn(\vo\,v) = \^N(yo,y^v)f 




cxp 



I/a 



(37) 



This probability distribution conserves probability along the classical trajectory y oc 
fk(v) when the condition (2iV + l)|F(/c)| ^> 1 is satisfied. As we shall show below in 
section 4, the requirement that the latter condition should be better and better satis- 
fied with the expansion of the Universe picks up a small subset of all possible classical 
trajectories fk{v) m phase space (namely, the quasi-isotropic mode of perturbations). 
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3 Heisenberg representation 



Let us perform now an independent analysis of the transition to semiclassical be- 
haviour in the Heisenberg representation. The crucial point here is that, as shown 
in 0, in the limit of infinite squeezing the mode functions fk(v) can be made real 
by a time- independent phase rotation (of course, in addition, |/&| — > oo). Then the 
momentum mode functions gk{v) become purely imaginary. As a result, the field 
operators y(k, 77) and p(k, rj) take the following form: 

y{k,r}) = f k (a(k,7 ?0 )+a t (-k, ?7o )) , (38) 
PiKv) = ~Wk (a(k,77o) +a t (-k,r/ )) = -i vfcv) • (39) 

So, in the limit |rj.| — > oo, the non-commutativity of y(k, rj) and p(k,rj) (and of all 
other operators, too) may be neglected. Due to this reason, we can introduce equiva- 
lent classical stochastic quantities yk(v) = fk e(k) where e(k) are c- number complex 
stochastic time-independent variables obeying the same statistics as y(k, r] ) at the 
time 770- They have zero average and unit dispersion: (e(k)) = 0, (e(k)e*(k / ) = 
^( 3 )(k — k'), though they are generally not Gaussian. The definition of the equiv- 
alence is the following: we say that a complex classical stochastic field 2/(77) with a 
corresponding momentum p(rj) and a probability distribution in phase space p(\y\, \p\) 
are equivalent to quantum operators 7/(77) an d P(v) if they give the same average values 
for an arbitrary function G(y,p): 

H (N\G(y,p)&(y,p)\N) H = J J dy 1 dy 2 dp 1 dp 2 p(\y\, \p\) \G(y,p)\ 2 (40) 

(here we omit k for brevity but we restore hats over operators). 

Of course, such a condition cannot be fulfilled generally, since quantum mechanics 
is not equivalent to classical field theory. However, in our case of infinite squeez- 
ing, due to relation (|39D between y and p (while the same relation between classical 
stochastic quantities y and p holds), eq. (flUD reduces to a much simpler form: 



H (N\G(y)&(y)\N) H = / / d yi dy 2 p(\y\) G(y) (41) 



where G(y) = G(y, p(y)) (since all operators are mutually commuting, we do not 
write any more hats over them). For the iV-particles initial state, the probability 
density function p(\y\) is given by 




P(\y\) = -ZH^Ll ( rjU ) exp I - Jfj- ) . (42) 
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So, in the limit of infinite squeezing and after use is made of relation ( p9[ ) the equiv- 
alence condition takes a form which is identical to the definition of an average value 
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of an operator in quantum mechanics with the probability distribution p equal to the 
modulus square of a wave function satisfying the Schrodinger equation. Thus, Q¥ip 
can be immediately proved by referring to the results of the previous section. How- 
ever, if we want to produce a completely independent derivation in the Heisenberg 
representation, we should not use the fact that p is determined by some solution of 
the Schrodinger equation but instead we have to make an actual calculation of both 
sides of eq. fl4T|) and prove their equality. 

Further, we omit tilda over G. We compute first the l.h.s. of fl4l|) : 



H (N\G(y)G\y)\N) H = £ q m q* n H (N\y m (k, V )y^(k, V )\N) H . (43) 

m,n=0 



We use now the expression (|39|) in the high squeezing limit when fk(v) can be made 
real by a time-independent phase rotation, and write: 

oo rm+n 

H (N\G(y)G\y)\N) H = £ q^j^ <0|a"(-k, Vo) a N (k, rj ) x 

m,n=0 V iV V H 

(a(k, 770) + at(-k, m )) m ( a t(k, rto) + a(-k, a tAr (k, r ?0 )a^(-k, ry )|0) H . (44) 
We find after a somehow tedious but straightforward calculation: 



H (N\G(y) G\y)\N) H = £ \ qm \ 2 f k m (v) £ '" H(A " 4 ^ 



m=0 i=sup(0,m-N) il ijn - i)\ (N - T71 + l)! 

(45) 

As for the r.h.s. of fl4i"D , we get: 

<Wlfep(lvl)|G(v)| a 
E // dyidy 2 -^ L N y m y* n ex P . (46) 



m,n=0 "1/* 

Replacing L^r by its explicit polynomial expression, we find: 

/ J d Vl dy 2 p (\y\ 2 ) \G(y)\ 2 

00 ^ (— 1V +J 'A7"! 2 /-oo 
- V I/7 l 2 f 2m V 1 ^ / - r m + i +ip- a; r /- r 

m =o i J=0 * ! J ! (iV-^)l(A^-j)! Jo 

^ , ,2 f2m f (-l)' +J >! 2 (m + »+j)! 
" io' gml/fc ^0 z! 2 i! 2 (iV-z)!(iV-i)! ' (4?) 

One finally checks, as we did numerically, that the sums in ( pE5|) and fl4T| ) are indeed 
equal: 

™ m! 2 (iV + *)! = » (-l)^iV! 2 (m + z+j)! 

i= ^-iV)^(^-i)! 2 (iV-m + i)! ^0 i! 2 j! 2 (iV-i)!(iV-i)! " 
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This ends the proof of equality (f41|). We note that the calculation done above yields 
in particular all the correlation functions. Indeed, if we set G(y) = y m , we get: 

«<A1|yP'"|iv>« = IAr t M m ^S!^ m + (V • < 49 > 

i=sup(0,m~N) % - \ m % Y- l iV m + l)\ 

We see also that all the following correlation functions do vanish: 

H (N\y m y^\N) H = if m^n. (50) 

Having in mind the comparison of the probability density obtained here with a Gaus- 
sian one, the two-point and four-point correlation functions are of particular impor- 
tance: 

H (N\\y\ 2 \N) H = (2N + l)\f k \\ 

H (N\\y\ 4 \N) H = (6N 2 + 6N + 2)\f k \ i . (51) 

The above results can be directly used in order to calculate the kurtosis Q (the 
irreducible 4-order moment) of the probability distribution: 

(N, ^^(k,^)! 4 ^, rj) If 1 \ 

V " (N, 77|Iz/(k, 77 )|2|iV, V y " 2^ (2N+iy) ■ 1 ) 

It vanishes, as should be the case, when N = and tends to —1/2 for iV — > oo. We 
will come back to this point in Section 6 and in the conclusion. For comparison, we 
present the corresponding results for a one-mode iV-particles state for which 

tp = (p(t)a + <p*(t)a\ [a,a f ] = l. (53) 

An example of such a state is the k = mode. Then 

= (2iV + l)|^| 2 , 
(if 4 ) = (6iV 2 + 6iV + 3)M 4 , 

° s w-^-^-ptfW)- (54) 

These expressions will be also useful for the calculation of the curtosis of y in the 
coordinate representation, see section 7. 



4 Conditions for semiclassical behaviour of pertur- 
bations 

Let us now analyze in more detail the conditions imposed on the mode functions fk{f]) 
of the quantum field y or on the equivalent stochastic functions yk(v) — /fc(^)e(k) by 
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the requirement of high squeezing \r k \ — > oo which is crucial for the quasi-classical 
equivalence as we have shown in sections 2 and 3. In the long-wavelenght limit 
k <C aH, the leading term in the mode functions fk(v) an d 9k(v) have the form [^, [?]: 

f k = d(k)a + C 2 (k)a f V JjL, g k = O^kA + ^1 (55) 

Joe a [Tj ) ±± a 

where Ci(k) is made real and positive by a phase rotation. We assume here that 
the integral converges at the lower limit that takes place for all physically interesting 
equations of state of matter at late times, in particular, if a(t) oc t n with n > 1/3). 
C*i(k) and C 2 (k) should satisfy the normalisation condition 

C^C 2 = - l - . (56) 

The first term in ([55]) is the quasi-isotropic mode, the second is the decaying mode. 
In the simplest inflationary models with one slow-rolling scalar field, a more concrete 
expression for C\{k) and C 2 (k) is obtained || [5J: 

Hi. ik 3/2 
Ci = -FV, C 2 = -^U- (57) 

where H k is the value of H{rj) at the moment of the first Hubble radius crossing during 
the inflationary stage. However, generally we need not restrict ourselves to this case 
only. The first term in the expression for g k in eq.(|55]) has a rather complicated 
integro-differential structure that can be defined by solving (0) using a perturbative 
expansion in powers of k 2 . The given estimate is valid for a wide class of smooth 
functions a (tj) including arbitrary power-law and quasi-de Sitter ones. 

Now let us insert ( j55|) into the definition of F(k), eq.(§). It is argued in p] that 
it is the condition |i^(fc)| ^> 1 that is required for the quasi-classical behaviour of 
the vacuum initial state. In the case of iV-particles state, the necessary condition is 
weaker, namely (2iV + l)|F(fc)| > 1 as shown above. Assuming a power- law or a 
quasi-de Sitter behaviour of the scale factor a (77) and omitting numerical coefficients 
of the order of unity, we obtain the following condition: 



If the solution is dominated by the quasi-isotropic mode, and the term with C 2 may 
be neglected, the condition ( |5*5D takes the form 

In terms of metric tensor perturbations (gravitational waves) which are related to 
(p = y/a by the relation 

h i:i = — ^ = V327TG <f etj (60) 
a 
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where ejj(k) is the polarisation tensor normalised by the condition e^e 1 - 7 = 1, i, j = 
1, 2, 3, we get: 

k 3 hl(k) ee (fc%(k)^*(k)> » ^ ■ (61) 

where lp = y/G is the Planck length and A = 2nak~ 1 is the wavelength of the 
perturbation. 

The condition fl5"9|) is increasingly better satisfied with the expansion of the Uni- 
verse. Thus, the quasi-isotropic mode becomes more and more classical at later times. 
For the simplest inflationary initial conditions (^7[), the inequalities (|59| ) and ( |6"TD are 
satisfied if XH(2N + 1) ^> 1, in particular, everywhere in the long-wavelenght regime 
XH > 1 even for iV = 0. 



Eqs. (59) and (61) were derived under the assumption XH 1. In the opposite 
case Aif -C 1, i.e. after the last Hubble crossing, eq. (54) of paper M must be used 
instead of (|53). This leads to the classicality condition: 



» | ■ ^r-r (62) 



which was derived in [|7j for the case N = 0. On the other hand, let us consider the 
case where the quasi-isotropic mode is negligible with respect to the decaying mode 
while the latter is imaginary ( 1 3^ C*2 1 "C |^yC 2 | in the representation where QCi = 0). 
This requires, in particular, C\ to be much less than the value given by fl5"7|) at the 
moment of any Hubble radius crossing. Note that is uniquely defined by C\ 
through the normalisation condition (|5(J) while 3ftC*2 is not restricted at all. Then the 
classicality condition takes the form: 

|C 2 | 2 » tftffQ » A L_ . ( 63 ) 

1 21 2N+V 9 X 3 H 2N + 1 y J 

These inequalities become weaker with the expansion of the Universe for all interesting 
cases (in particular, for a(t) oc t n with n > |) and may even cease to be valid. Thus, 
the imaginary decaying mode becomes less classical and more quantum at later times. 
Finally, a sufficiently large real and positive decaying mode (3?C 2 ^> \QC 2 \ = -^r) is 
always classical, since in this case the inequality (^) is already satisfied due to the 
third term in the left-hand side. In contrast, a negative decaying mode (KC2 < 0, 
1 3F£C 2 1 ^> ISC2I) has a rather peculiar behaviour: being classical at early and late 
times, it ceases to be classical at the moment when the first and the third term in 
the left-hand side of ( |58|) almost cancel each other. 

All that was said above refered either to the case of a massless minimally cou- 
pled scalar field (p, or to tensor metric perturbations (gravitational waves). In the 
case of scalar (or adiabatic) metric perturbations, one should first either specify a 
gauge, or choose some gauge-independent quantity describing the perturbations. Let 
us consider the gauge-independent quantity $ which coincides with the Newtonian 
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potential in the longitudinal gauge (g 0Q = 1 + 2$). However, one should keep in mind 
that semiclassical behaviour may be more pronouced in terms of some quantities and 
less pronouced, or even absent, for some other quantities. 

Let us further restrict ourselves to the case where pressure perturbations are diag- 
onal. This includes all cases where matter is a mixture of an arbirary number of scalar 
fields with arbitrary potentials (but minimally coupled to the gravitational field) and 
an arbitrary number of ideal fluids. Then, as is well known, in the long-wave regime 
A A j ~ H~ x cs where Aj is the Jeans lenght and cs is the sound velocity, a Fourier 
mode of the gravitational potential may be written as: 

$ fc = C\(k) (l - - [\ 2 d V ) - ^f-C 2 (k)- (64) 
\ a Jrn J At a 

where C*i(k) can be made real and positive by a phase rotation (see, e.g., [15| for 



a simple derivation of a more general statement valid in the synchronous gauge). 
Once more, the first term in the right-hand side of ( |64| ) is the quasi-isotropic mode 
(which may be really called the growing mode now since the corresponding density 
perturbation increases with the expansion of the Universe), the second term is the 
decaying mode, rji should be as small as possible, but not less than the moment of 
the first Hubble radius crossing at the inflationary stage. 
If the operator $(k, if) is represented as: 



Phi(k,rj) = <f> k (rj) a(k,77o) + ®* k (v) a f (-k, 



then the normalisation condition for the coefficients Ci(k) and C 2 (k) of the mode 
functions &k{v) i n the long- wave regime has the same form as 



C&C 2 = -~ . (66) 

For the case of matter in a Universe consisting of one minimally coupled scalar field 
(p, this statement follows from the fact that in this case, canonical variables analogous 



to y and p are given by [|TB| : 



Z = a( = a5<p L + ^ , = (67) 

ti a 

where 6^p L is the scalar field perturbation in the longitudinal gauge. Using two exact 
relations between and $&: 

AtxGo 2 (p 1 (H V 

1 H ( a V . . 

& = sc^(s*0 (68) 
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and the exact background equation H = — 47rG0 2 , it is straightforward to obtain the 
leading terms in ( k and in the long-wavelenght region from eq. 



Cl H + ° 2 H 



1 R2 J 

oo (f z a d 



H 



H 



o ^ 



C?77 . 



(69) 



Since, on the other hand, should satisfy the normalisation condition — 



z, we arrive to (|66|). Note also the exact equal-time commutation relation following 



from the canonical commutation relations between £ and ir^ and eq. (|6~8): 



$(k,77),$ t (k,77) 



.4vrG # 



fc 2 a 



(70) 



For the simplest inflationary models driven by one slow-rolling minimally coupled 
scalar field <p: 



C 1 (k) 



V2k?/*\<p k \ 



C 2 (k) 



lk^ k \ 



(71) 



where H k and ip k are the values of H(t) and (p(t) at the moment of the first Hub- 
ble radius crossing during the inflationnary stage. Eq. (|7l| ) leads to the standard 
expressions for density perturbations generated during inflation. 

Let us now turn to the case of matter consisting of one hydrodynamical component 
with an ideal fluid equation of state P = P(e), where P is the pressure, e is the energy 
density. The canonical variables analogous to y and p in this case were first found by 
Lukash [|i~7f : 



q = aq 



Tin 



3a H 2 ( a x 
— $ 



V8nG aH \H 



c s V 3 ' Cs 



dP 

de ' 



, a' _ 
q - —q 
a 



(72) 



(it is assumed here that c$ ^ 0). If the sound velocity c$ is constant, then a oc a. 
The corresponding wave equation for phonon mode functions has the form: 



a" 



(73) 



Note also the second relation between q k and 



2 — Ha q k 



k 2 



(74) 



17 



and two useful identities for the background metric in this case: 

2H o H 



1 _L r 2 



(75) 



The corresponding leading long-wavelength behaviour of is: 



Qk 



Ci 



87rG ~ /■* (it c ; 



Co ^r- 



1 + f 



(76) 



Using (|72|) and (ff5|), it is easy to verify that the commutation relation (|70 ) holds in 
this case too, which leads to (|66"D once more. Though we derived the commutation 
relation (|70|) for two particular cases (one scalar field or one ideal fluid), we believe 
that it has a much wider range of validity and we leave the proof of this assertion for 
another publication. 

The condition for semiclassical behaviour of the gravitational potential is: 



(^(M) $ f (k,??) + c.c.) > [$(k,r/), $ + (k,7/)] 



(77) 



Substituting here the long-wave asymptote ( |64|) and assuming hydrodynamical matter 
with a constant sound velocity (that leads to the power-law behaviour a(t) oc a n with 
| < ft < |), we get the following condition of classicality similar to fl58|) : 



in 



4ttG off 3 1 1 v A; 2 a 

If the growing mode is dominating, then 

~ 2 4vrG a/7 3 1 
Ci > 



> 



1 



2N + 1 



(7* 



P Pc 2 o 2iV + 1 ' 



(79) 



This condition is increasingly better satisfied with the expansion of the Universe. In 
terms of the gravitational potential, this leads to 



k^\k) = (fc 3 |$(k)| 2 ) > 



ll\H ? - 



2N 



10) 



On the other hand, for a purely decaying imaginary mode (IKC2I <C ISC2I), we get 



\C 2 2 » 



4ttG H 2N + 1 



F$ 2 (A;) > 



Z 2 # 



A 2iV + l 



This condition weakens with the expansion of the Universe. Note that the conditions 
for the gravitational potential ( [SOUSl]) are more restrictive than the corresponding 
conditions for gravitational waves ( pl|j63| ) since the right-hand sides of the former 
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conditions are in (XH) 2 ^> 1 times larger (however, the growing mode is always semi- 
classical in the long-wavelength regime for the most interesting case of perturbations 
produced during inflation with the initial conditions (|7l|)). As pointed above, this 
peculiar feature may be a result of $ not being the best quantity to study quantum- 
to-classical transition. 

In conclusion, we see that the condition for a perturbation to become more and 
more classical at later times (i.e. with the expansion of the Universe) chooses a 
subclass of all solutions for perturbations, namely, the quasi-isotropic modes. We will 
return to this point at the end of the paper. 



5 The Wigner function 

The Wigner function was proposed as a candidate for the probability density function 
of a quantum mechanical system in phase space. According to the rules of quantum 
mechanics such a probability density can not exist, so one can only hope to come as 
close as possible to it by requiring from good candidates, like the Wigner function, 
to fulfill some basic properties of the probability density, but not all of them. In 
particular, the Wigner function is real, but not definite positive. Let us calculate the 
Wigner function for the system under consideration : 

W NN = W (y ,yo,Po,P* ) 

= T^ja / / fci dx 2 e-^ + ^V* N (y - * N (y + | jf7 ) .(82) 

We will use the simplified notations y(k, rjo) = y — y\ + iy 2 and the same convention 
for p(k, ?7o). The quantity p\ , resp. p 2 , is canonically conjugate to y% , resp. y 2 , and so 
pi = 2 pi, p 2 = 2 p 2 . 

We consider first the Wigner function W a @ of the coherent initial state ^ a p. Ac- 
tually, it is clear from ( |17D that it will be similar to the Wigner function Wqq of the 
vacuum |0, rj). The following result is obtained 



7T 



M s L.m \ s L m y \ (83) 



I x |2 a t \ I £ 

\h\ ) V \Ji 



k\ 



We recognize the vacuum Wigner function Wqq with p — > p — (p) and y —>■ y — (y). 

We now turn our attention to the Wigner function W^n of A^-particles states. 
Our method will be to compute the generating functional of W^n- F° r this purpose 
we first write down the generating functional of W mn : 



J] W mn z m v n = — j j dx 1 dx 2 e~^ + ^ x 



m,n=0 
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m,n=0 

The generating functional of the Laguerre polynomials is given by: 

00 1 rz 

Y: L m (x)z m = exp(--f-) (\z\ < 1) . (85) 

m=0 



We insert the expression for the vacuum wave function into fl8"3|), replace each sum 
of Laguerre polynomials by the generating function ( j35|) and rearrange terms in the 
exponentials. The integration reduces to Fourier transforms of gaussians. In order to 
give a compact result, we will set: 

The result of the integration can then be written in the following way: 

^ rrr m n 11 ( (l^| 2 + l^| 2 )(l + ^)\ 

X, W mn z m v n = — exp' 1 ' 



m,n=0 



7T 2 1 — ZV \ 1 — ZV J 

Y\ 2 + \P\ 2 + 21^^! + 2zy 2 P 2 )t/ 



exp 



exp 



1 — ZV 

\Y\ 2 + \P\ 2 - 2iY 1 P 1 - 2iY 2 P 2 )z" 
1 — zv 



11/ (\Y\ 2 + \P\ 2 )(l + zv) 



■ exp 



7T 2 1 — ZV 1 V 1 — ZV 

00 , iy.2 + ip|2 + 2 ^y 1 p 1 + 2iY 2 P 2 ) i v i 



X 



i=o 

~ (-|F| 2 + \P\ 2 - 2iY 1 P 1 - 1iY 2 P 2 )i Z i 

§ • (87) 

We have calculated the finite form of the series J2mn=o W mn z m v n , but what we really 
want in order to get W^n is the finite form of the series J2m=o ^mm,u m - This means 
that we must select terms of equal order in z and v in the double sum and replace the 
product zv by u. Looking at the last expression fl37|), we see that the terms that must 
be kept are those with i — j in the two sums. We can perform the transformation 
and write the relevant generating functional: 



m=0 1X 2 1-U \ 1-U 



X 



~ ((|F| 2 - |P| 2 ) 2 + (2Y 1 P 1 + 2Y 2 P 2 ) 2 ) l u i 
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7T 1 — U \ 1 — U I 



Ji 



u ((\Y\ 2 - |P| 2 ) 2 + (2Y 1 P 1 + 2Y 2 P, 



2 ) 2 



1/2 ' 







1 — M 



where Jo is the Bessel function. 

It is now possible to identify the generating functional of a product of Laguerre 



polynomials by using the identity (eq. 8.976.1 in [|Tj]) 



m=0 

If we take: 



\TTf\Tf^m 1 ( (z + y)A j ( 2(-xyu) 1/2 ) 

2^ L m (x)L m (y)u =- exp J — . (89) 

1 — it \ 1 — u I \ 1 — u 



x = \Y\ 2 + \P\ 2 + 2Y 1 P 2 - 2Y 2 P 1 = \Y - iP\ 2 , 

y = \Y\ 2 + \P\ 2 - 2Y 1 P 2 + 2Y 2 P 1 = \Y + iP\ 2 , (90) 



we can finally write ( jS8|) as: 

OO 1 oo 

£ W mmu m = -e-d y l 2+ l p l 2 ) J2 L m {x)L m {y)u m . (91) 

m=0 71 m=0 

This gives the final result: 

W mm = ±e~^ 2+ ^ L m (\Y -iP\ 2 ) L m (\Y + iP\ 2 ) 

7T Z 

= W m L m {\Y-iP\ 2 )L m {\Y + iP\ 2 ) . (92) 

Transition to semiclassical behaviour can be analyzed using this result. We see from 
the last expression that the Wigner function, again in the limit of infinite squeezing 
\tk\ — > oo, will be delta-concentrated around classical trajectories of the system. 
Indeed, V^oo; the Wigner function corresponding to the initial vacuum satisfies |7j 

W 00 - l^ol 2 5 (P) (93) 

and, therefore, also 

W NN -> |^v| 2 5 (P) (94) 

where S(P) is the 5-distribution. Thus, we obtain one more independent derivation of 
the relation (^) between y and p in the limit \rk\ — * oo which results in the reduction 
of © to ©. 

An interesting point is that, in contrast to the vacuum initial state where the 
Wigner function is positive definite, here the Wigner function is no longer positive 
definite before the limiting transition |rfc| — > oo leading to (0) is performed. Thus, 
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W cannot be interpreted as a probability distribution in classical phase space for 
any finite \rk\- This property actually constitutes the core of the criticism of the 
Wigner function approach in the papers |TT]. However, we do not want to insist on 



this classical interpretation of the Wigner function before the limiting transition (why 
should we?), and we use the Wigner function only to show that, for high squeezing, it 
becomes concentrated near a bunch of classical trajectories in phase space with initial 
stochastic amplitude y(k,r] ) obeying the probability density function (|37|). For this 
reason, too, we deliberately use the word "concentrated" instead of "peaked" because 
the latter term, though often used may create a wrong impression about the positivity 
of the fine structure of the Wigner function around these classical trajectories. On 
the other hand, after the limiting transition \rk\ — > oo the Wigner function acquires 
the structure (|94]) which may be used as a classical probability distribution in phase 
space. Moreover, this probability distribution is sufficient for all problems where 
oscillating quantum fine structure ("wave pattern") around classical trajectories can 
be neglected. 

This discussion also sheds some new light on the question as to what extent any 
coarse graining is necessary to obtain the transition to semiclassical behaviour of 
quantum cosmological perturbations. A coarse graining becomes really necessary if 
one intends to obtain some classical probability distribution in phase space from the 
Wigner function before the limiting transition \r^\ — > oo. However, for quantities for 
which the above mentioned quantum fine structure is unimportant (and all quantities 
observable in forseeable future belong to this class) all parameters resulting from 
the correct (i.e. physically justified) coarse graining disappear after the transition 
\rk\ — > oo, so we are allowed to forget about coarse graining at all. As a result, for 
these quantities, decoherence (understood as a partial loss of quantum coherence, 
with the remaining part of coherence being expressable in classical terms) is reached 
without consideration of any actual mechanism of coarse graining and decoherence 
(that is why we call it "decoherence without decoherence"). On the other hand, there 
do exist quantities and problems for which this fine structure is important (e.g., it 
is essential for the calculation of the entropy of cosmological perturbations). But 
then one should not introduce a coarse graining simply by hand or average over some 
degrees of freedom of the Universe arbitrarily called as "unobservable" . All such 
choices have to be justified by consideration of real physical processes. 

The very specific case of a coherent initial state, for which the Wigner function is 
given by eq. (|83|) , is still instructive since it clearly shows a smooth transition between 
deterministic and purely stochastic behaviour of perturbations in the semiclassical 
limit. Indeed, if |(y(^))| 2 3> max-fi, | /&(??) | 2 }, we get an effective deterministic 
classical behaviour (a time-dependent Bose condensate of the field y) . Realization of 
this inequality at all times is possible, e.g., by choosing initial conditions (y) and (p) 
proportional to fk(vo) an d 9k(Vo) resp. with a large coefficient of proportionality. In 
the opposite limit \fk(v)\ 2 ^ max{| (y{r))) | 2 , ^} characteristic for real cosmological 
perturbations, the average value of perturbations is negligible with respect to their 
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rms value. Then quantum perturbations become equivalent (in the sense defined 
above) to classical quantities with a stochastic amplitude, a fixed phase defined by 
the 5 multipliers in the expression for the Wigner function ([83]), and a practically zero 
average value. Of course, generally both average and rms values may be non-zero. 
We see that the appearance of deterministic quantities in the semiclassical limit is a 
very specific phenomenon restricted to a narrow class of both special initial states and 
special time evolution (with not too many "particles" created). In particular, even the 
choice of a coherent initial state which is already classical and deterministic at the 
initial moment 7]o does not guarantee that it remains so at later times. A constructive 
counterexample is provided by the case where the initial conditions (y)o and (p)o 
correspond to the pure decaying mode in the long-wavelenghth regime XH ^> 1. On 
the contrary, transition to effective stochastic classical behaviour is a generic effect 
which arises as a result of high squeezing of a very wide class of initial quantum states. 



6 Thermal state 

Until now we have considered N-particles initial states. Let us now consider the 
important case of a thermal initial state. Such a state arises for example in the 



cosmological scenario where the stage of "thermal inflation" [|18j takes place after 
relaxation of the radiation-dominated Universe to the state of thermal equilibrium if 
we mean by i]q the moment when this stage of late inflation begins. In this case, the 
probability density in the y(k)-representation is given by : 



7>(y(k),y(-k)) = (y\(Y,e~ N/No \N)(N\)\y) 

\N=0 I 

oo 

= A(N ) Ee"^;(!/(k))fe(!/(k)) . (95) 

where Nq is the ratio of the temperature to the energy field quanta for given k and 
A(Nq) is a normalization constant. So, actually Nq depends on k but we drop the k 
for brevity of notation. The probability to find N quanta in the initial state drops 
rapidly for N > N . 

This probability density can be calculated using the generating functionnal of 
squared Laguerre polynomials, already given in 

V{\y\)) = ^f^H te-^'Ll 

Ljv=o 






y\ 


2 2e~ 1 ' A 


1/* 


2 i _ e -l/N 




X I j 1 2 eX P 



k 
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= 7Tw exp {-& coth ^ Jo (&^k) ' (96) 

When iVo <S 1, i-e. the temperature is much smaller than the energy of the quanta 
for given k , we see that V is almost Gaussian, since the exponential dominates the 
Bessel function. This is of course to be expected since then our thermal state almost 
reduces to the vacuum state, there are practically no thermal excitations. What can 
be said if N grows, from what temperature on will we get a noticeable departure 
from a Gaussian probability density ? It turns out that for any value of N , the 
probability V is a symmetric single-peaked function centered around zero, at first 
sight quite similar to a gaussian. Therefore, the relevant discriminating parameter 
to be calculated in order to study the gaussianity of this distribution is the kurtosis. 
Fortunately, the latter can be found analytically. Since V depends only on \y\, its 
even momenta are given by the following integrals: 







y\ 2n V(\y\)27r\y\d\y\ 



roo 

\f k \ 2n x"exp(-xcoth — -) J 

JO Z1\q 



IX 
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Four different integrals of the Bessel functions are given in [14|, 6.621.1, with the 



results expressed in terms of the hypergeometric function F. Using the third one, we 
check the normalisation Iq = 1. With the second one, we calculate I2 and I4, and 
finally the kurtosis Q: 



h H/*| 2 coth^F(4o;l;cosh-^)H/*| 2 coth^, (9* 
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Q = h /(I 2 y - 2 = cosh' 2 — - . (100) 

Looking at fig.[l], we see that there is a sharp transition between a low-temperature 
gaussian behaviour (Q = 0) and a high-temperature non-gaussian behaviour (Q — 1). 
Intermediate values of the kurtosis imply that N is close to one, within only one order 
of magnitude. Fig.^| shows how the gaussian transforms into another distribution at 
high temperature. Scales of cosmological interest were inside, most of them even deep 
inside, the Hubble radius at the onset of inflation. Their energy is therefore much 
higher than the temperature at that time. This corresponds to values of N Q <C 1 for 
which, as can be seen from ( 100 ) and from fig.^], the probability looks very much like 
a Gaussian. 
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7 Conclusions and discussion 



We have shown here that, for a wide class of initial non-vacuum states including 
all iV-particles, thermal and most of the coherent states, quantum inhomogeneous 
cosmological perturbations (both adiabatic perturbations and gravitational waves) 
generically acquire semiclassical behaviour in the regime when their wavelenghth is 
much larger than the Hubble radius. In the language of quantum field theory, the 
latter regime corresponds to the high squeezing limit |rjt| — > oo. Then an equivalent 
classical description of these perturbations in terms of c-number stochastic quantities 
becomes possible. This, according to our paradigm, constitutes the transition from 
quantum to effectively classical behaviour of cosmological perturbations. Thus, we 
generalize the results of obtained for the perturbations of the vacuum initial state. 
The main difference between vacuum and non-vacuum initial conditions is that the 
statistics of the stochastic amplitude of the equivalent classical field modes is Gaussian 
in the former case but generally non-Gaussian in the latter case. 

We proved this equivalence in three different independent ways. First, we ex- 
plicitly constructed the wave function in the Schrodinger representation thereby ex- 
hibiting its WKB shape. Then, we have shown that the operators in the Heisenberg 
representation correspond to classical stochastic functions in that limit. Finally, we 
have computed the Wigner function and shown that it becomes concentrated along 
classical trajectories. Since we do not in general choose a unique vacuum state from 
an infinite number of unitary non-equivalent vacuum states at the initial moment 
r] (the case of the inflationary scenario where a unique choice is possible was used 
as an illustration only), our results are valid for iV-particles, thermal and coherent 
states constructed using any possible vacuum state. From a mathematical point of 
view, this ambiguity is reflected in the freedom to choose the constant parameters 
Ci(k) and C^k) characterizing the long- wavelength regime (|53D which satisfy only 
one normalization condition fl51|). 

It turns out that the generalization from vacuum to non-vacuum initial conditions 
does not produce any new problem of principle for the quantum-to-classical tran- 
sition. The only thing which is important for it is that the classicality conditions 
for perturbations derived in section 4 should be satisfied (in particular, it is enough 
to have a sufficienly large amplitude of the quasi-isotropic mode of perturbations in 
the long- wavelength regime). The Wigner function for iV-particles or thermal ini- 
tial states (but not for coherent initial states) is not positive definite for a finite r&. 
However, we have shown that this does not prevent the transition to the semiclassical 
behaviour in the limit |r^| — > oo expressed by eq. (P4|). We believe that this addresses 
the criticism against the transition to the semiclassical behaviour raised in [|TTJ since 
the generalization of our considerations to the case where the scale factor a(rj) is 
itself quantized but still in the WKB regime does not bring new principally different 
features. 

Note also that we never use the so called "semiclassical gravity" , or one-loop ap- 
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proximation, = 8^(2(2"^) to study the evolution of quantum perturbations and 
the quantum-to-classical transition for them (these equations are usually understood 
as equations for the average value (g^) of the metric tensor). This approximation 
is inappropriate for the problem under consideration, not only because one-loop ap- 
proximation is not exact when only a few quantum fields are important, but mainly 
due to the fact that these equations are grossly insufficient in the typical case arising 
in cosmology where average values of perturbations (Sgik) are negligible (or even ex- 
actly zero) as compared to their rms values, and all their higher moments should be 
calculated (see also the discussion at the end of section 5). So, one should go beyond 
the one-loop, or the Gaussian, approximation. 

Now we discuss the question whether this non-Gaussianity of perturbaions which 
appears in the case of non-vacuum initial states is observable, and if it can be used 
to exclude such states using observational data on ^f- angular fluctuations of the 
CMB temperature. First, one should take into account that observable quantities in 
physical space are real and consist of sums of quantities of the type (y(k, r])e* kr + c.c), 
irrespective of the fact whether y is an operator or a c-number stochastic quantity. 
As a result, the statistics of quantities in real space is the same as for one- mode 
states. In particular, in the Gaussian case the probability distribution is the one- 
dimensional Gaussian distribution, not the two-dimensional one that prevails for the 
complex amplitude y(k, if). Then the value of the curtosis Q for quantities in real 
space in the case of an iV-particles initial state is given by the expression (|54|)(e.g., 

ip can be taken as 4p itself). Previously, the skewness S = ((^) ) an d topological 
characteristics, like the genus, of the CMB temperature fluctuations were mainly 
considered. However, S = for all iV-particles initial states. So, we propose to 
measure the kurtosis Q since this is more interesting. For non-vacuum states which 
we considered, Q is negative. 

For a non-Gaussian distribution, it is easy to see that Q lies below the cosmic 
variance level of purely Gaussian perturbations for all initial non-vacuuum states 
where different k modes are ^-correlated and rms values of perturbations depend on 
|k| only. This results due to the averaging over angular directions which is equivalent 
to summing over all k modes with the same |k|. Such a sum will rapidly become 
Gaussian due to the central limit theorem. In particular, the curtosis of the sum of 
n equally distributed modes with kurtosis Q is equal to Q/n ~ 1/n while the cosmic 
variance limit for the Gaussian distribution ~ \j\fn. Still, the non-zero kurtosis 
can be observable for less symmetric and more correlated initial states. An extreme 
example is when AT f^'^ can take only 2 values ±Ao with equal probability. Then 
the kurtosis Q = —2 and it certainly lies above the cosmic variance level. On the 
other hand, this discussion shows that there is a wide class of initial non-vacuum 
states which cannot be distinguished from the vacuum initial state by looking at the 
statistics of observational quantities, so it is not easy to disprove them. 

Finally, we have shown the crucial role played by the quasi-isotropic mode of 
perturbations in the quantum-to-classical transition. Namely, if this mode is present 
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then, as time goes on and the Universe expands, perturbations become more and 
more semiclassical. On the contrary, this does not take place for the decaying mode. 
Note that there is no quasi-isotropic mode in the case of the exactly homogeneous 
and isotropic FRW model (even fully quantized) or of the cosmological model of the 
Bianchi I type. There are known problems with the quantum-to-classical transtion 
in FRW quantum cosmology (even for large values of the scale factor) summarized 



in paper |19fl . We will not discuss here to what extent these problems are really 
dangerous for the FRW cosmology, but we would like to point out that the results of 
our paper suggest that the final solution to these problems will require consideration of 
more general, inhomogeneous cosmological models which contain the quasi-isotropic 
mode of perturbations (or its generalization to the fully non-linear case). Also, it is 
tempting to connect the arrow of time in the Universe with the quantum-to-classical 
transition. Then a possible relation between the arrow of time and the growth of 
inhomogeneities in the Universe (the transition from order to chaos) emerges. 
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A 



Figure 1: The kurtosis Q as a function of log(iVo) for a squeezed thermal state specified 
by N (for given k). 
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Figure 2: The probability distribution as a function of log(iV ) again for a squeezed 
thermal state. In order to visualize the shape of the distribution and not its width, 
we plot V, for given N , as a function of \y\/ \J < y y^ >. We see that for N <C 1, it 
is basically Gaussian. 
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